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OLEG KARPENKOV 

Abstract. Recently we showed that Hessenberg matrices are proper to represent con- 
jugacy classes in SX(n, Z). In this paper we focus on the reducibility properties in the 
set of Hessenberg matrices of SL(3, Z). We investigate the first interesting open case 
here: the case of matrices having one real and two complex conjugate eigenvalues. 
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Introduction 

In current paper we study the SX(3,Z) integer conjugacy classes via reduced represen- 
tatives. Recall that two matrices Mi and M2 in SL(3, Z) are integer conjugate if there 
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exists a matrix X in GL(3, Z) such that 

M 2 = .Y.U|.Y ; . 

For matrices with a pair of complex conjugate eigenvalues we discovered the following 
phenomenon: Hessenberg matrices distinguish corresponding conjugacy classes asymptoti- 
cally (Theorem 2.6). We show that a similar statement is not true for the case of operators 
with three real eigenvalues. 

Background. In classical approach to SL(n, Z)-conjugacy problem one splits SL(n, Q)- 
conjugacy classes into SL(n, Z) conjugacy classes. After that the problem is reduced to 
certain problems related to orders of algebraic fields extended by the roots of characteristic 
polynomial of the corresponding matrices (like computing their class numbers, etc.). 

In [18] we introduced multidimensional analog of Gauss Reduction Theorem, it is an 
alternative approach to the conjugacy problem (for two-dimensional Gauss Reduction 
Theory we refer to [17], [26], and [27]). In multidimensional Gauss Reduction Theory the 
key role play Hessenberg matrices, they generalize reduced matrices in Gauss Reduction 
Theory. Hessenberg matrices are matrices that vanish below the superdiagonal (for more 
information see in [31]). They appear in the work [13] by K. Hessenberg for the first 
time, they were later used in QR- algorithms ([12], [32], [30]). In [18] we defined a natural 
notion of Hessenberg complexity for Hessenberg matrices, which is a nonnegative integer 
function, and showed that each integer conjugacy class of irreducible matrices has only 
finite number of Hessenberg matrices with minimal Hessenberg complexity. 

Description of the paper. We start in Section 1 with general definitions and notation. 
Further in Section 2 we formulate main results of current paper: Theorem 2.1 on parabolic 
structure of the sets of Hessenberg matrices with two complex conjugate eigenvalues, and 
Theorem 2.6 on asymptotic reducibility of matrices in these sets. In Section 3 we prove 
Theorem 2.1. Further in Section 4 we show some necessary tools that we use in the proof of 
Theorem 2.6 (Markoff-Davenport characteristic, Klein- Voronoi continued fractions, etc.). 
Then in Section 5 we give a proof of Theorem 2.6. Finally in Section 6 we formulate 
several open problems. 

Acknowledgment. The work is partially supported by FWF grant M 1273-N18. The 
author is grateful to E. I. Pavlovskaya and H. W. Lenstra for useful remarks. 

1. Definitions and notation 

1.1. ^-reduced matrices. A matrix M of the form 

an 0-12 ai3 
a 2 \ a 2 2 0,23 
a 32 a33 

is called an (upper) Hessenberg matrix. We say that the Hessenberg type of M is 

(an, fl2i|ai2, 022, ^32). 
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Definition 1.1. A Hessenberg matrix in SL(3, Z) is said to be perfect if we have 

< an < a 21 \ 
< a 12 < a 32 ; 
< a 22 < a 32 . 

Definition 1.2. The Hessenberg complexity of a Hessenberg matrix M of Hessenberg type 
(aii,a 2 i\ai 2 ,a 22 ,a,3 2 ) is the number ai 2 a 23 , we denote it by 

Notice that ?(M) equals the volume of the parallelepiped spanned by the following 
vectors v = (1,0,0), M(v), and M 2 (v). 

Definition 1.3. We say that a perfect Hessenberg matrix M is ^-reduced if its Hessenberg 
complexity is the least possible. Otherwise we say that the matrix is q-nonreduced. 

In [18] we proved the following result. 

Theorem 1.4. i). Any conjugacy class o/S'L(3,Z) contains a ^-reduced matrix, 
ii). The number of q -reduced matrices is finite in any integer conjugacy class. 

The results of current paper give evidences concerning the fact that the majority of 
Hessenberg matrices with two complex conjugate and one real eigenvalues are ^-reduced. 

1.2. Perfect Hessenberg matrices of a given Hessenberg type. In this paper we 
study three-dimensional perfect Hessenberg SX(3, Z)-matrices with irreducible character- 
istic polynomials. There are two main geometrically essentially different cases of SX(3, Z)- 
matrices: the real spectrum (or RS- for short) case when the characteristic polynomials of 
matrices have only real eigenvalues, and the nonreal spectrum (or NRS- for short) case of 
matrices with a pair of complex conjugate and one real eigenvalues. 

Denote the set of all «SX(3, Z)-matrices of a fixed a Hessenberg type Q by H(Q). Let 
H V (Q) be the subset of all NRS-matrices in NRS(fl). 

Definition 1.5. Let Q = (an, a 21 \a 12 , a 22 , a 32 ). Consider v = (ai 3 , a 2 3, a 33 ) such that the 
determinant of the matrix (a^) equals 1. Denote 

(an a 12 a u m + a 12 n + a 13 
a 2i a 22 a 21 m + a 22 n + a 23 
a 32 a 32 n + a 33 

It is clear that 

H(Q) = {H£(m,n)\m eZ,neZ} 

Here to choose v means to choose the origin O in the plane H(Q). So the set H(Q) 
has the structure of two-dimensional plane. We denote by OMN the coordinate system 
corresponding to the parameters (m,n). 

Let T>^(m,n) denote the discriminant of the characteristic polynomial of H^(m,n). 
Then the set NRS(fl) is defined by the following inequality in variables n and m: 

V v n (m,n)<0. 
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in = —n 
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x(-i) = o 



rn 



m = n — 2 



x(i) = o 



Figure 1. The family of matrices of Hessenberg type (0, 1|0, 0, 1). 



Example 1.6. In Figure 1 we show the subset of NRS-matrices NRS({0, 1 10, 0, 1)). For 
this example we choose v = (0, 0, 1). 



We start in Subsection 2.1 with the formulation of a supplementary theorem on para- 
bolic structure of the set of NRS-matrices, we give the proof later in Section 3. Further 
in Subsection 2.2 we formulate the main result on asymptotic uniqueness of ^-reduced 
matrices, the proof is shown in Section 5. In Subsection 2.3 with describe examples of 
families of matrices with fixed Hessenberg type. 

2.1. Parabolic structure of the set of NRS-matrices. The set NRS({0, 1|0,0, 1)) on 
Figure 1 "reminds" the set of points with integer coordinates in the union of the convex 
hulls of two parabolas. Let us formalize this in a general statement. 

Consider the matrix Hq(0,0) = (a^) and define b±, b 2 , and 6 3 as coefficients of charac- 
teristic polynomial of this matrix in variable t: 



In the case of SX(3,Z) we have 6 3 = 1, nevertheless we write 63 for generality reasons. 
For the family H^(m,n) we define the following two quadratic functions 



2. Formulation of main results and examples 



-t 3 + M 2 - ht + 63. 



Pi,n(m, n) = m - a^ 2 - fan - 71; 



n /a 2 im — a u n 

p 2 ,n{m, n) = a 2 [ 

a 2 i V a 21 




a 2 \m — dun 



) 



-72, 
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where 



7i 



Q32 

4a 2 i 
an - a 2 2 



^33 



2a 2i 
4fr 2 - g 

4a 2i a 32 



a 2 

ft 
72 : 



032&21 

463 

o 2 



26, 



4&1&3 



4a 2 ia 32 & 3 



Denote by Br{0) the interior of the circle of radius R centered at the origin (0, 0) in 
the real plane OMN of the family H^(m,n). For a real number t we denote 

A t = {(m,n) I (p lt n(m,n) - t)(p 2 ,n(m, n) -t) < 0}. 

Theorem 2.1. For any positive e there exists R > such that in the complement to 
B R {0) the following inclusions hold 

A £ c NRS(n) c A_ £ . 
We give a proof of this theorem in Section 3. 

2.2. Theorem on asymptotic uniqueness of ^-reduced NRS-matrices. A point is 
called integer if all its coordinates are integers. A ray is said to be integer if its vertex is 
integer and it contains integer points distinct to the vertex. 

Definition 2.2. An integer ray in H(Q) is said to be an NRS-ray if all its integer points 
correspond to NRS-matrices. A direction is said to be asymptotic for the set NRS(Q) if 
there exists an NRS-ray with this direction. 

As it is stated in Theorem 2.1, for any Hessenberg type Q the set NRS(Q) almost 
coincides with the union of the convex hulls of two parabolas. This implies the following 
statement. 

Proposition 2.3. Let fl = (a u , a 2 i|ai 2 , a 22 , a 32 ). There are exactly two asymptotic direc- 
tions for the set NRS(Q), they are defined by the vectors (—1,0) and (ou,o 2 i). □ 

Let us consider a Hessenberg type Q = (a n , a 2 i\a 12 , a 22 , a 32 ) and an appropriate integer 
vector v . 

Definition 2.4. Consider a family of Hessenberg matrices Hq. Denote 



RZn, v = {H v n(rn-t,n)\teZ> o y, 

R T,a,v = {H&(m+a n t,n+a 2 it)\t G Z> }. 



By R™Q V (t) or respectively by R^Q V (t) we denote the t-th element in the corresponding 
family. 

Remark 2.5. The families R™q v and R™q v coincide with the sets of all integer points of 
some rays with directions (—1,0) and (an,a 2 i) respectively. Conversely, from Proposi- 
tion 2.3 it follows that the set of integer points of any NRS-ray coincides either with R™q v 
or with R™n v for some integers m and n. 
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Figure 2. Any NRS-ray contains finitely many ^-nonreduced matrices. 

On Figure 2 we show in dark gray two NRS-rays: Ri^f 2 \i 13) (00-1) f rom the an d 
^<i72|i,i,3),(o,o,-i) from the ri S ht - 

Now we are ready to formulate the main result on asymptotic behavior of NRS-matrices, 
we prove it later in Section 5. 

Theorem 2.6. (On asymptotic ^-reducibility and uniqueness.) i). Any NRS-ray 
(as on Figure 2) contains only finitely many q-nonreduced matrices, 
ii). Any NRS-ray contains only finitely many matrices that have more than one integer 
conjugate ^-reduced matrix. 

Example 2.7. Any NRS-ray for the Hessenberg type (0, 1|0, 0, 1) contains only ^-reduced 
perfect matrices. Experiments show that any NRS-ray for (0, 1 1 1, 0, 2) contains at most 
one ^-nonreduced matrix (see in Figure 3 on page 7). 

2.3. Examples of NRS-matrices for a given Hessenberg type. In this subsection 
we study several examples of families NRS(Q) for the Hessenberg types: 

(0,1|0,0,1), (0,1|1,0,2), (0, 1|1, 1,2>, and (1,2|1,1,3). 

In Figures 3, 4, and 5 the dark gray squares correspond to ^-nonreduced matrices. We 
also fill with gray the squares corresponding to ^-reduced Hessenberg matrices that are 
n-th powers (n > 2) of some integer matrices. 

Hessenberg perfect NRS-matrices H^'^ ^(m, n) . The Hessenberg complexity of 
all these matrices is 1, and, therefore, they are all ^-reduced, see the family on Figure 1 
on page 4. 

Hessenberg perfect NRS-matrices H^'^ 2 Jm, n) . The Hessenberg complexity of 
these matrices equals 2. Experiments show that 12 of such matrices are ^-nonreduced, 
see the family in Figure 3. It is conjectured that all others Hessenberg matrices of 
^^((0,111,0,2)) are ^-reduced. 
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m+n = — 1 




Figure 3. The family of Hessenberg matrices H^'^ 02 ^(m,n). 



m+n = — 1 




Figure 4. The family of Hessenberg matrices \ 2 ^ (m, n) 



Hessenberg perfect NRS-matrices H^'^\ 2 ^(m,n). The Hessenberg complexity of 
these matrices equals 2. We have found 12 ^-nonreduced matrices in the family. It is 
conjectured that all other Hessenberg matrices of NRS({0, 1|1, 1,2)) are ^-reduced. See 
in Figure 4. 

Hessenberg perfect NRS-matrices H^!^ 3 Jm, n). This is a more complicated exam- 
ple of a family of Hessenberg perfect NRS-matrices, their complexity equals 12. We have 
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x(-i) = o 




Figure 5. The family of Hessenberg matrices H^'^i 13 y 

found 27 ^-nonreduced matrices in the family. It is conjectured that all other Hessenberg 
matrices of NRS({1,2\1, 1,3)) are ^-reduced. See in Figure 5. 

3. Proof of Theorem 2.1 
We start the proof with several lemmas, but first let us give a small remark. 
Remark. The set NRS(Q) is defined by the inequality 

V v n (m,n)<0. 

In the left part of the inequality there is a polynomial of degree 4 in variables m and 
n. Note that the product lQa 21 a\ 2 b^ (pi,n(m, n)p2,n(m, n)) is a good approximation to 
T>Q(m, n) at infinity: the polynomial 

V^(m,n) - lQal 1 al 2 b 3 (p hn (m,n)p 2tn (m,n)) 

is a polynomial of degree 2 in variables m and n. 

Lemma 3.1. The curve 2^o'i|oo i)( m > n ) = ^ s contained in the domain defined by the 
inequalities: 

( (m 2 -4n + 3)(n 2 + 4m + 3) > 

\ (m 2 -An- 3)(n 2 + 4m - 3) - 72 < 

Remark. Lemma 3.1 implies that the curve f[o'i|oo i>( m ' n ) = is contained in some 
tubular neighborhood of the curve 

(m 2 - 4n)(n 2 + 4m) = 0. 
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Proof. Note that 

^(0,110,0,1) 1 



f/o'iioo i)( m ' n ) = i 171 " 2 ~~ 4^)(n 2 + 4m) — 2mn — 27. 



Thus, we have 

p <Moki)( m ' n ) ~ ~ An + 3 )( n2 + 4m + 3 ) = -2(n - 3) 2 - 2(m + 3) 2 - (n + mf < 0, 



and 

X){^[° ) 01) (m,n)-(m 2 -4n-3)(n 2 + 4m-3) + 72 = 2(n-3) 2 + 2(m + 3) 2 + (n-m) 2 > 0. 

Therefore, the curve f |o'i|oo i)( m ' n ) = is contained in the domain defined in the lemma. 

□ 

Lemma 3.2. For any Q = (an, a2i|ai2, Q22, 032) ^ere exists an affine {not necessarily 
integer) transformation of the plane OMN taking the curve V^(m,n) = to the curve 

Proof. Let Hq(0,0) = (a i: j). Note that a matrix H^(m,n) is rational conjugate to the 
matrix 

by the matrix 

(1 an a n + ai 2 a 2 i 
a 2 i 011021 + o 2 ia 22 
a 2 ia 32 

Therefore, the curve V^m, n) = is mapped to the curve X^o'ifo i)( m > n ) = bijectively. 
In OMN coordinates this map corresponds to the following affine transformation 

m \ ^ / a 2i a 32 m - ana 32 n \ / a 23 a 32 - ana 33 + a i2 a 2 i - a 22 a 33 - ana 22 
n / \ 032^ / \ an + ^22 + a 33 

This completes the proof of the lemma. □ 

Proof of Theorem 2.1. Consider a family of matrices pi^O, i) + e,t) with real 

parameter t. Direct calculations show that for e 7^ the discriminant of the matrices for 
this family is a polynomial of the forth degree in variable t, and 

V v n (- Pl ,n(0,t) +e,t)= \a 21 a\ 2 et A + 0(f). 

Therefore, there exists a neighborhood of infinity with respect to the variable t such that 
the function £>^(— pi ^(0, t) + e,t) is positive for positive e in this neighborhood, and 
negative for negative e. 

Hence for a given e there exists a sufficiently large iV\ = iV\ (e) such that for any t > Ni 
there exists a solution of the equation V^(m,n) = at the segment with endpoints 

(-Pi,a(0,t)+e,t) and ( - pi,n(0, t) - e, t) 

of the plane OMN. 
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Now we examine the family in variable t for the second parabola: 
(t-a llP2 , a (t, 0)- °" 2 e, -a 21 p 2 , n (t, 0)~^= 

V v a ii +a 2i v a ii +a 



By the same reasons, for a given £ there exists a sufficiently large N 2 = N 2 (e) such that 
for any t > N 2 there exists a solution of the equation D^(m,n) = at the segment with 
endpoints 

{t-a llP2 , n (t, 0)-^m=e, -a 21 p 2 , n (t, 0)-^m=e) and 
of the plane OMN. 

We have shown that for any of the four branches two parabolas defined by pi t n(m, n) = 
and p 2t n{ m , n ) — there exists (at least) one branch of P^(m,n) = contained in the 
e-tube of the chosen parabolic branch if we are far enough from the origin. 

From Lemma 3.1 we know that 2^o'i|oo i)( m ' n ) = is contained in some tubular 
neighborhood of 

Pl,(0,l|0,0,l)( m > n )P2,(0,l|0,0,l>(^,^) = 0. 

Then by Lemma 3.2 the curve V^(m, n) = is contained in some tubular neighborhood 
of the curve 

Pi,n(m,n)p2,n(rn,n) = 

outside some ball centered at the origin. Finally, by Viet Theorem, the intersection of the 
curve T>n(m, n) = with each of the parallel lines 

an + Q21 , 

t t : n — m = t 

a 2 i 

contains at most 4 points. Therefore, there exists sufficiently large T such that for any t > 
T the intersection of the curve T>^(m, n) = and i t contains exactly 4 points corresponding 
to the branches of the parabolas pi,n(m, n) — and p 2j n(m, n) = lying in A_ e \ A e . 

Hence, there exists R = R(e, N\, N 2) T) such that in the complement to the ball B R (0) 
we have 

A e C NRS{Q) C A_ £ . 

The proof of Theorem 2.1 is completed. □ 

4. Supplementary tools for the proof of Theorem 2.6 

In this section we introduce several notions that we use in the proof of Theorem 2.6. In 
Subsection 4.1 we introduce Markoff-Davenport characteristic that represents the Hessen- 
berg complexity. Further in Subsection 4.2 we show how to construct perfect Hessenberg 
matrices (M\v) conjugate to a given one. Finally in Subsection 4.3 we give the definition of 
Klein- Voronoi continued fractions, formulate a theorem on construction of ^-reduced oper- 
ators via vertices of a fundamental domain of the corresponding Klein- Voronoi continued 
fraction, and prove one supplementary statement on geometry of continued fractions. 
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4.1. MD-characteristics. The study of the Markoff-Davenport characteristics is closely 
related to the theory of minima of absolute values of homogeneous forms with integer 
coefficients in n-variables of degree n. One of the first works in this area was written 
by A. Markoff [28] for the decomposable forms (into the product of real linear forms) for 
n — 2. Further, H. Davenport in series of works [7], [8], [9], [10], and [11] made first steps 
for the case of decomposable forms for n = 3. 

Consider A G SX(n, Z). Denote by P(A,v) the parallelepiped spanned by vectors v, 
A(v), A n ~\v), i.e., 

f n—l 

P(A,v) = lO + Y,*iA i (v) 

^ z=0 

where O is the origin. 

Definition 4.1. The Markoff-Davenport characteristic (or MD- characteristic, for short) 
of an SL(n, Z)-operator A is a functional: 

A A : R n ->■ R defined by A A (v) = V(P(A, v)), 

where V(P(A,v)) is the nonoriented volume of P(A,v). 

Remark 4.2. Consider an operator A with Hessenberg matrix M in some integer basis. 
Then the Hessenberg complexity q(M) equals the value of MD-characteristic A j4 (l,0,0). 

We continue with the following general definition. 

Definition 4.3. The group of all GL(3, Z)-operators commuting with A is called the 
Dirichlet group and denoted by S(A). 

For MD-characteristic we have the following invariance property. 

Proposition 4.4. Consider A e SL(n, Z) and let B e 2(A). Then for an arbitrary v we 
have 

A A (v) = A A (B(v)). 

Basically, this means that the MD-characteristic naturally defines a function over the 
set of all orbits of the Dirichlet group. 

4.2. Construction of a perfect Hessenberg matrix (M\v) conjugate to a given 

one. Let us show how to construct perfect Hessenberg matrices integer conjugate to a 
given one. 

Algorithm to construct perfect Hessenberg matrices. 

Input Data. We are given by an SL(3, Z)-matrix M of an operator A with irreducible 
characteristic polynomial over Q and an integer primitive (i.e., with relatively prime 
coordinates) vector v. 

Step 1. We put t\ = v. 



< A; < l,i = 0, . ..,n-l 
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Step 2. Choose an integer primitive vector of the plane spanned by v and A(v) on the 
minimal possible nonzero Euclidean distance from the line spanned by v, denote it by g 2 . 
Find the coordinates qn and a 2 i from the vector decomposition 

M e l) = <?ll e l + °2102- 

Find 6n and a n as integer quotients and reminders: 

qn = \a 2 i\bu + a n . 

Define 

e 2 = sign(a 2 i)fi'2 + h^. 

Step 3. Choose an integer primitive vector g% e M 3 on minimal possible nonzero Eu- 
clidean distance from the plane spanned by e 1 and e 2 . Find the coordinates q± 2 , q 22 , and 
a 32 from the vector decomposition 

A(e 2 ) = q\ 2 e x + g 22 e 2 + a 32 g3- 

Find 612, b 22 , a i2 , and a 22 as integer quotients and reminders: 

qi2 = |a 32 |6i 2 + a 12 and g 22 = \a 32 \b 22 + a 12 . 

Then we have 

e 3 = h 2 ei + b 22 e 2 + sign(a 32 )g 3 . 

Output Data. Let C be a transition matrix to the basis {e 1 ,e 2 ,e 3 }. In the output we 
have the perfect Hessenberg matrix CMC^ 1 . 

Definition 4.5. Consider an SL(3, Z)-matrix M with irreducible characteristic polyno- 
mial over Q and an integer primitive vector v. Starting from M and v the above algorithm 
generates a perfect Hessenberg matrix, we denote it by (M\v) 

Remark 4.6. In [18] we showed that any perfect Hessenberg matrix integer conjugate to 
M is represented as (M\v) for a certain integer primitive vector v. 

4.3. Klein- Voronoi continued fractions. In the proof of Theorem 2.6 we essentially 
use the geometric construction of Klein- Voronoi continued fractions. In [19] and [20] 
F. Klein proposed a multidimensional generalization of continued fractions to totally real 
case. First attempts to find analogous construction in other cases were made by G. Voronoi 
in his dissertation [33]. In 1985 J. A. Buchmann in his papers [5] and [6] proposed to use 
Voronoi's extension to compute of fundamental units in orders. We use a slightly modified 
definition of Klein- Voronoi continued fraction from the paper [18]. 

4.3.1. RS-case. Let us first briefly recall Klein's definition of two-dimensional continued 
fraction in totally real case. Consider an operator A in GL(3, Z) with three real distinct 
eigenvalues. This operator has three distinct invariant planes passing through the origin. 
The complement to the union of these planes consists of 8 open orthants. Let us choose 
an arbitrary orthant. 
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Definition 4.7. The boundary of the convex hull of all integer points except the origin 
in the closure of the orthant is called the sail. The set of all 8 sails of the space 1R 3 is 
called the 2 -dimensional continued fraction in the sense of Klein. 

For further information on Klein continued fractions we refer to the following litera- 
ture: [24], [3], [1], [2], [21] [23], [22], [25], [14], [15], [16], [29] etc. 

4.3.2. NRS-case. Consider an operator A in GL (3, M) with distinct eigenvalues. Suppose 
that it has a real eigenvalue r and complex conjugate eigenvalues c and c. 

Denote by T X (A) the set of all real operators commuting with A such that they have 
a real eigenvalue equals 1 are with absolute value of both complex eigenvalues equal one. 
Actually, T l (A) is an abelian group with operation of matrix multiplication isomorphic 
toS 1 . ' 

For v G M 3 we denote 

T A (v) = {B(v)\BeT 1 (A)}. 

If v is a real eigenvector then Ta(v) consists of one point. Otherwise (in general case) 
T A {v) is homeomorphic to S 1 . 

Let g\ be a real eigenvector with eigenvalue r; g 2 and g 3 be vectors corresponding to 
the real and imaginary parts of some complex eigenvector with eigenvalue c. Consider 
the coordinate system OXYZ corresponding to the basis Denote by 7r the (/tri- 

dimensional plane OXY, and by n + — the half-plane of n defined by y > 0. 

Proposition 4.8. For any v the orbit T A (v) intersects the half-plane ir + in a unique 
point. □ 

Definition 4.9. A point p G ir+ is said to be ^-integer if the orbit T A (p) contains at least 
one integer point. 

The invariant hyperplane x = of operator A divides 7r + into two arcwise connected 
components. 

Definition 4.10. The convex hull of all tt- integer points except the origin contained in 
a given arcwise connected component is called a factor-sail of the operator A. The set of 
both factor-sails is said to be the factor- continued fraction for the operator A. 
The union of all orbits T A (*) in R n represented by the points in the factor-sail is called 
the sail of the operator A. The set of all sails is said to be the continued fraction for the 
operator A in the sense of Klein- Voronoi (see in Figure 6 below). 

It is clear that the factor-sail is a broken line. The corresponding sail is the surface of 
elliptic rotation of the factor-sail around the eigenline of A. The cones corresponding to 
rotation of edges (vertices) are called factor-edges (factor-vertices) . 

4.3.3. Algebraic continued fractions. Consider an operator A in GL(3, Z) with a real eigen- 
value r and two complex conjugate distinct eigenvalues c and c. Suppose also that the 
characteristic polynomial of A is irreducible over Q. 
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The Dirichlet group (of GL(3, Z)-operators commuting with A) takes the Klein- 

Voronoi continued fraction to itself but maybe exchange the sails. By Dirichlet unit 
theorem (see in [4]) the Dirichlet group S(A) is always homomorphic to Z © Z/2Z. 

Definition 4.11. A fundamental domain of the Klein-Voronoi continued fraction for A 
is a collection of open orbit-vertices and orbit-edges such that for any orbit-face F of the 
continued fraction there exists a unique orbit-face F' in this collection and an operator 
T e E(A) such that F = T(F'). 

Example 4.12. Consider an operator 



It has one real and two complex conjugate eigenvalues. In Figure 6a we show in light gray 
the halfplane n + , the invariant plane for A corresponding to complex conjugate eigenvalues 
is colored in dark gray. The boundary line of the halfplane ir + is an invariant line of A, it 
contains real eigenvectors of A. 

The halfplane tt + is shown in Figure 6b. The invariant plane intersects 7T + in a ray 
separating n + into two connected components. A point of n + is colored in black if and 
only if it is a 7r-integer point. The boundaries of the convex hulls in each part of 7r + arc 
two factor-sails. Notice that, one factor-sail is taken to another by the induced action of 
the operator —Id, where Id is the identity operator. 

In Figure 6c we show one of the sails of Klein-Voronoi continued fraction for A. There 
are three visible orbit-vertices, they correspond to integer vectors (1,0,0), (0,1,0), and 
(0, 0, 1): the large dark points (0, 1, 0) and (0, 0, 1) are visible on the corresponding orbit- 
vertices, the point (1,0,0) is on the backside of the continued fraction. 

A fundamental domain of the operator consists of one orbit- vertex and one orbit edge. 
For instance, one can take the orbit-vertex corresponding to the point (1, 0, 0) and the 
orbit-edge corresponding to the "tube" connecting orbit-vertices for the points (1,0,0) 
and (0,1,0). 

In the proof of Theorem 2.6 we use the following result on construction of <;-reduced 
operators via vertices of fundamental domains of Klein-Voronoi continued fractions. 

Theorem 4.13. ([18]) Consider an SX(n, Z) = operator A with matrix M having distinct 
eigenvalues. Let U be a fundamental domain of the Klein- Voronoi continued fractions for 
A. Then we have: 

(i) For any ^-reduced matrix M integer conjugate to M there exists v e U such that 
M = (M\v). 

(ii) Let v e U. The matrix (M\v) is ^-reduced if and only if the MD- characteristic 
Aa(v) attains its minimal value. □ 

4.3.4. One general fact on fundamental domain of Klein-Voronoi continued fractions for 
NRS-matrices of <SX(3, Z). Further we use the following statement. 
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a) 



b) 



Figure 6. A Klein- Voronoi continued fraction: a) the cone 7r + and the 
eigenplane; b) the continued factor-fraction; c) one of the sails. 



Consider an NRS-operator A in SL(3, Z) and any integer point x distinct from the 
origin. Denote by T° A (p) the convex hull of the union of two orbits corresponding to the 
points p and A(p). For any integer k we denote by T k A (p) the set A k (T° A (#)). 

Proposition 4.14. Let A be an NRS-operator in SX(3,Z) and p be an integer point 
distinct from the origin. Then there exists a fundamental domain of the Klein- Voronoi 
continued fraction for A with all (integer) orbit-vertices contained in the set T° A (p). 

The proof is based on the following lemma. Let 

rA( P ) = \Jr k A ( P ). 



Lemma 4.15. Consider A e SX(3,Z) with NRS-matrix and let p be any integer point 
distinct from the origin. Then one of the Klein- Voronoi sails for A is contained in the set 

Proof. Notice that the set T A (p) is a union of orbits. Let us project T A (p) to the halfplane 
7r + . The set T A (p) projects to the closure of the complement of the convex hull for the 
points n(A k (p)) for all integer number k in the angle defined by eigenspaces. Since all 
the points A k (p) are integer, their convex hull is contained in the convex hull of all points 
corresponding to integer orbits in the angle. Hence tt(Ta(p)) contains the projection of 
the sail. Therefore, the set T A (p) contains one of the sails. □ 

Proof of Proposition 4-H- Since —Id exchange the sails, one can choose a fundamental 
domain entirely contained in one sail. Let this sail contains a point p. By Lemma 4.15 
Ta{p) contains this sail. Therefore all the orb it- vertices of a fundamental domain for 
Klein- Voronoi continued fraction can be chosen from the factor-set T A (p). □ 
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5. Proof of Theorem 2.6 

5.1. Geometry of Klein- Voronoi continued fractions for matrices of R™q v . Let 
us show the following statement. 

Proposition 5.1. Consider an NRS-ray R™q v - Then there exists C > such that for 
any t > C there exists a fundamental domain for the Klein- Voronoi continued fraction 
of the matrix R™n v (t) such that all integer points in this domain are contained in the 
triangle with vertices (1,0,0), (an,a 2 i,0), and (-an,-a 2 i,0). 

We begin with the case of matrices of Hessenberg type Qq — (0, 1|0, 0, 1). Such matrices 
form a family H(Qq) with real parameters m and n as before: 

^(o,'i|'o,o,i)( m ' ri ) = 
Here v = (1,0,0). 

Lemma 5.2. Let R™q v be an NRS-ray. Then for any e > there exists C > such 
that for any t > C the convex hull of the union of two orbit-vertices 

T^ o(t) (l,0,0) and ^,^(0,1,0) 

is contained in the e-tubular neighborhood of the convex hull of three points (1,0,0), 
(0,1,0), (0,-1,0). 

Remark. Actually Lemma 5.2 means that the corresponding domain tends to be flat while 
the parameter t tends to infinity. 

Proof. Let us find the asymptotics of eigenvectors and eigenplanes for operators -R™^ „ (t) 
while t tends to +oo. Denote the real eigenvector of R™q vo (t) by e(t). We have 

e (t)=^((i,o,o) + o(r 1 )) 

for some nonzero real /i. 

Consider the unique invariant real plane of the operator R™j£ VQ (t) (it corresponds to 
the pair of complex conjugate eigenvalues). Notice that this plane is a union of all orbits 
T R m,n ( t )(w) for arbitrary vectors w of this plane. Any such orbit is an ellipse with axes 

Xg m3/X (t) and Xg m m(t) for some positive real number A, where 

= (0,t,0) + O(l), 
^min(0 = (0,0,t 1 / 2 )+O(r 1 /2). 

Actually, the vectors 5 f m ax(^)±^fi , min(^) are two complex eigenvectors of R™n Vo (t) . For 
the ratio of the lengths of maximal and minimal axes of any orbit we have the following 
asymptotic estimate: 

X\9mUt)\ = | t |l/2 + 0( | t |-l/2 )- 
A|0minW| 
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Since 

(1,0,0) --e(t) = 0(|t|- 1 ), 
fx 

the minimal axis of the orbit-vertex T R ™,n ,,Jl 0,0) is asymptotically not greater than 

0(t~ r ). Therefore, the length of the maximal axis is asymptotically not greater than 
some function of type 0(\t\~ 1 / 2 ). Hence, the orbit of the point (1, 0, 0) is contained in the 
(Ci |t| — 1 / 2 )-ball of the point (1,0,0), where C\ is a constant that does not depend on t. 
We have 

(0,1,0) -^ max (t) = 0(|t|- 1 ). 

Therefore, the length of the maximal axis of the orbit-vertex T R ™,g ( t )(l, 0, 0) is asymp- 
totically not greater than some function 1 + 0(t~ l l 2 ). Hence, the length of the minimal 
axis is asymptotically not greater than some function 0(|t| -1 / 2 ). This implies that the 
orbit of the point (0, 1, 0) is contained in the (C2|£| _1//2 )-tubular neighborhood of the seg- 
ment with vertices (0,1,0) and (0,-1,0), where C 2 is a constant that does not depend 
on t. 

Therefore, the convex hull of the union of two orbit- vertices 

T *DUM (1 '°' 0) and T KS (hV( M°^ 0) 

is contained in the C-tubular neighborhood of the triangle with vertices (1, 0, 0), (0, 1, 0), 
(0, —1, 0), where C = max(Ci, C 2 )|t|~ 1 / 2 . This concludes the proof of the lemma. □ 

Let us now formulate a similar statement for the general case of Hessenberg matrices. 

Corollary 5.3. Let Q = (an, a2i|ai2, G22, 032) and R™q v an NRS-ray. Then for any 
e > there exists C > such that for any t > C the convex hull of the union of two 
orbit-vertices 

7 ^fr,(*)( 1 ' ' ) and T R^(t)( a H' a 2i>o) 

is contained in the e-tubular neighborhood of the convex hull of three points (1,0,0), 
(an, O2i,0), (-an, -a 2 i,0). 

Proof. Denote Q = (an, a 2 i|ai 2 , a 22 , 032) and choose 

( a 21 a 32 — a 32 a ll a ll a 22 ~ a 21 a 12 
a 32 -an - a 2 2 

1 

A direct calculation shows that 

^(0,'iioai) (^i(m,n) - ~^ ^ 2 ( n °)) = XH ^( m ~ t,n)X~ 1 , 

where l\ and l 2 are linear functions with coefficients depending only on an, 021, 012, 
a 2 2, and a 32 . Therefore, the ray R™q v after the described change of coordinates and a 

homothety is taken to the ray R™q Vq of matrices with Hessenberg type f2 = (0, 1|0, 0, 1) 
for certain rh and n. 
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Lemma 5.2 implies the following. For any e > there exists a positive constant such 
that for any t greater than this constant the convex hull of the union of two orbit- vertices 

T r?£ (*)(!> °>°) and T R m,n ( ) (0,1,0) 

is contained in the e-tubular neighborhood of the triangle with vertices (1, 0, 0), (0, 1, 0), 
(0,-1,0). 

Now if we reformulate the last statement for the family of operators in old coordinates, 
then we get the statement of the corollary. □ 

Proof of Proposition 5.1. We note that the operator R™Q v (t) takes the point (1,0,0) 
to the point (an, a 2 i, 0). Therefore, the convex hull of the union of two or bit- vert ices 

^Dv^ 1 ' ' ) and T «r^w( aii ' a2i '°) 

(we denote it by W(t)) coincides with the set T° R m,n ^(1,0,0). 

From Proposition 4.14 it follows that there exists a fundamental domain for the con- 
tinued fraction with all its orb it- vert ices contained in W(t). Choose a sufficiently small 
e such that the e -tubular neighborhood of the triangle with vertices 

(1,0,0), (01,1,02,1,0), and (-an, -a 2 i, 0) 

does not contain integer points distinct from the points of the triangle. From Corol- 
lary 5.3 it follows that for a sufficiently large t the set W(t) is contained in the eo-tubular 
neighborhood of the triangle. This implies the statement of Proposition 5.1. □ 

5.2. Geometry of Klein- Voronoi continued fractions for matrices of R™q v . Now 
let us study the remaining case of the rays of matrices with asymptotic direction (an, a 2 i). 

Proposition 5.4. Consider an NRS-ray R™q v - Then there exists C > such that for 
any t > C there exists a fundamental domain for the Klein- Voronoi continued fraction 
of the matrix R™n >v (t) such that all integer points in this domain are contained in the 
triangle with vertices (1,0,0), (—1,0,0), and (an,a 2 i,0). 

The proof of this proposition is based on the corollary of the following lemma. We 
remind that Q = (0, 1|0, 0, 1). 

Lemma 5.5. Let R™$i oVo be a> n NRS-ray. Then for any e > there exists C > such 
that for any t > C the convex hull of the union of two orbit-vertices 

T R ^ ( t)(l,0,0) and T R ^ (t) (0,l,0) 

is contained in the e-tubular neighborhood of the convex hull of three points (1,0,0), 
(-1,0,0), (0,1,0). 

Proof. First, we note that the continued fractions for the operators A and A^ 1 coincide. 
Secondly, the following holds: 

^(oJJai)^' n + t )= XH (o'mo,i)(- n - *' ~ m ) x ~\ 



ON ASYMPTOTIC REDUCTIBILITY IN SL{3, Z) 



19 



where 

/ -1 -n-t 
X= -1 -m 
\ -1 

Thus, in the new coordinates we obtain the equivalent statement for the ray RiQ~^(t). 
Now Lemma 5.5 follows directly from Lemma 5.2. □ 

Corollary 5.6. Let Q = (an, 021(012, 022, 032) and R™q v be an NRS-ray. Then for any 
e > there exists C > suc/i £/ia£ /or any t > C the convex hull of the union of two 
orbit-vertices 

(1,0,0) and T H v {m+ailttn+oait) (au,a 2 i,0) 

is contained in the e-tubular neighborhood of the triangle with vertices (1, 0, 0), (—1, 0, 0), 
and (an, 021, 0). □ 

Remark. We omit the proofs of Corollary 5.6 and Proposition 5.4, since they repeat the 
proofs of Corollary 5.3 and Proposition 5.1. 

5.3. Conclusion of the proof. Let us finally conclude the proof of Theorem 2.6. Let A 
be an operator with Hessenberg matrix M in SX(3, Z). By Theorem 4.13 any ^-reduced 
matrix congruent to M is constructed as the matrix (M\v), where v is an integer vector 
in an arbitrary chosen fundamental domain of the Klein- Voronoi continued fractions for 
A, in addition v should be the minimum of the absolute value of MD-characteristic on the 
integer lattice except the origin. To calculate ^-reduce matrices we find all such minima 
of MD-characteristics in appropriate fundamental domains. 

The case of NRS-rays with asymptotic direction (—1,0). Consider an NRS-ray R™£ v - 
By Proposition 5.1 there exists C > such that for any integer t > C we can choose a 
fundamental domain for the Klein- Voronoi continued fraction of R™n v (t) such that all its 
integer points are in the triangle with vertices (1, 0, 0), (an, 021, 0), and (—an, —021, 0). 

This triangle contains only finitely many integer points, all of them have the last coor- 
dinate equal to zero. The value of the MD-characteristic for a point (x, y, 0) equals: 

(a 2 \x - a xx y)a\ 2 yH + C, 

where the constant C does not depend on t, it depends only on x, y, and Q. Therefore, 
for any point (x, y, 0) the MD-characteristic is linear with respect to the parameter t, and 
it increases with growth of t. The only exceptions are the points of type A(1,0, 0) and 
°2i 5 0) (for integers A and //). The values of MD-characteristic are constant in these 
points with respect to the parameter t. 

Since there are finitely many integer points in the triangle (1,0,0), (an,a 2 i,0), and 
(— an, — a 2 i, 0), for sufficiently large t the MD-characteristic at points of the triangle 
attains the minima only at (1,0,0) and at (an,a 2 i,0). Since R™n v (t) takes the point 
(1,0,0) to the point (an, 021, 0), a fundamental domain may contain only one of these 
two points, let it be (1,0,0). 
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Therefore, for sufficiently large t the minimum of MD-characteristic at the integer points 
of the chosen fundamental domain is unique and it is attained at point (1, 0, 0). Hence by 
Theorem 4.13 for sufficiently large t the matrix 

#, (1 '°' 0) . Am-t.n) (1,0,0)) =H { , im , Jm-t,n) 

(an, a2i|ai2, 022,032) v 1 ) \ 1 1 ) J (an, a 2 i|ai2, 022,032) v ' > 

is the only ^-reduced matrix in the conjugacy class. This implies both statements of 
Theorem 2.6 for the ray -R^fiV 

Therefore, Theorem 2.6 holds for any NRS-ray with asymptotic direction (—1,0). 

The case of NRS-rays with asymptotic direction (an, 021). This case is similar to the case 
of NRS-rays with asymptotic direction (—1, 0), so we omit the proof here. 

Proof of Theorem 2.6 is completed. □ 



6. Open problems 

In this section we formulate open questions on the structure of the sets of NRS-matrices 
and briefly describe the situation for RS-matrices. 

NRS-matrices. As we have shown in Theorem 2.6 the number of ^-nonreduced matrices 
in NRS-rays is always finite. Here we conjecture a stronger statement. 

Conjecture 1. Let f2 be an arbitrary Hessenberg type. All but a finite number of NRS- 
matrices of type Q are ^-reduced. 

If the answer to this conjecture is positive we immediately have the following general 
question. 

Problem 2. Study the asymptotics of the number of ^-nonreduced NRS-matrices with 
respect to the growth of Hessenberg complexity. 



Denote the conjectured number of ^-nonreduced NRS-matrices of Hessenberg type f2 by 
Numerous calculations give rise to the following table for all types with Hessenberg 
complexity less than 5. 





(0,10,0,1) 


(0,11,0,2) 


(0,1 1,1,2) 


(0,1 1,0,3) 


(0,1 1,1,3) 


(0,11,2,3) 




1 


2 


2 


3 


3 


3 







12 


12 


6 


10 


10 




(0,12,0,3) 


(0,12,1,3) 


(0,12,2,3) 


(1,2 0,0,1) 


(0,11,0,4) 


(0,11,1,4) 




3 


3 


3 


4 


4 


4 




14 


10 


10 


94 


6 


8 




(0,1 1,2,4) 


(0,1 1,3,4) 


(0,13,0,4) 


(0,13,1,4) 


(0,13,2,4) 


(0,13,3,4) 




4 


4 


4 


4 


4 


4 




10 


8 


10 


12 


8 


8 



ON ASYMPTOTIC REDUCIBILITY IN SL{3, Z) 



21 



RS-matrices. We conclude this paper with a few words about real spectra matrices (i.e., 
about SX(3, Z) -matrices with three distinct real roots). Mostly we consider the family 
H((0, 1 1 1, 0, 2)), the situation with the other Hessenberg types is similar. 
Recall that 

/ 1 ra + 1 

ff (i'i|w>( m . fl ) =10 m 
\ 2 2n + 1 

This matrix is of Hessenberg type (0, 1|1, 0, 2), its Hessenberg complexity equals 2. Hence 
^(o il'i o 2)( m > n ) * s ^-reduced if and only if it is not integer conjugate to some matrix of 
unit Hessenberg complexity, all such matrices are of Hessenberg type (0, 1|0, 0, 1). 
In Figure 7 we show all matrices H^b^ Q2 Am,n) with 

-20 <m,n< 20. 

The square in the intersection of the m-th column with the n-th row corresponds to the 
matrix H^'^ 02 Jm,n). It is colored in black if the characteristic polynomial has rational 
roots. The square is colored in gray if the characteristic polynomial is irreducible and 
there exists an integer vector (x, y, z) with the coordinates satisfying 

-1000 < x,y,z < 1000, 

such that the MD-characteristic of ff( ,i|i,o,2>(l) 0, l)(m, n) equals 1 at (x,y,z). All the 
rest squares are white. 

If a square is gray, then the corresponding matrix is ^-nonreduced, see Remark 4.2. If 
a square is white, then we cannot conclude whether the matrix is ^-reduced or not (since 
the integer vector with unit MD-characteristic may have coordinates with absolute values 
greater than 1000). 

It is most probable that white squares in Figure 7 represent ^-reduced matrices. We 
have checked explicitly all the squares with 

-10 < m,n < 10. 

These matrices are contained inside the big black square shown on the figure. All white 
squares inside it correspond to ^-reduced matrices. 

We show a boundary broken line between the NRS- and .RS'-squares in gray. 

Remark. In Figure 7 the NRS-domain is easily visualized, it almost completely consists of 
white squares. While RS-domain contains relatively large number of black squares. This 
indicates a significant difference between RS- and NRS-cases. 

Direct calculations of the corresponding MD-characteristic show that the following 
proposition holds. 

Proposition 6.1. If an integer m+n is odd, then ^o'i|'io 2}( m > n ) ^ s ^-reduced. □ 

Remark. From one hand Proposition 6.1 implies the existence of rays entirely consisting 
of ^-reduced matrices. From the other hand in contrast to the NRS-matrices this is not 
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always the case for RS-matrices. For instance, all matrices corresponding to integer points 
of the lines 

l)m = n; 2)m = n + 2; 3)m = — n; 4)m = —n — 2; 5)n = 3m — 4; 6)m = 3n + 6 

are ^-reduced (we do note state that the list of such lines is complete). 

So Theorem 2.6 does not have a direct generalization to the RS-case and we end up 
with the following problem. 

Problem 3. What is the percentage of ^-reduced matrices among matrices of a given 
Hessenberg type f2? 

It is more likely that almost all Hessenberg matrices are ^-reduced (except for some 
measure zero subset). 

References 

[1] V. I. Arnold. Higher-dimensional continued fractions. Regul. Chaotic Dyn., 3(3):10-17, 1998. J. Moser 
at 70 (Russian). 

[2] V. I. Arnold. Preface. In Pseudo-periodic topology, volume 197 of Amer. Math. Soc. Transl. Ser. 2, 
pages ix-xii. Amer. Math. Soc, Providence, RI, 1999. 



ON ASYMPTOTIC REDUCTIBILITY IN SL{3, Z) 



23 



[3] V. I. Arnold. Continued fractions. Moscow: Moscow Center of Continuous Mathematical Education, 
2002. 

[4] A. I. Borevich and I. R. Shafarevich. Number theory. Translated from the Russian by Newcomb 
Greenleaf. Pure and Applied Mathematics, Vol. 20. Academic Press, New York, 1966. 

[5] J. Buchmann. A generalization of Voronoi's unit algorithm. I. J. Number Theory, 20(2):177-191, 
1985. 

[6] J. Buchmann. A generalization of Voronoi's unit algorithm. II. J. Number Theory, 20(2):192-209, 
1985. 

[7] H. Davenport. On the product of three homogeneous linear forms. I. Proc. London Math. Soc, 
13:139-145, 1938. 

[8] H. Davenport. On the product of three homogeneous linear forms. II. Proc. London Math. Soc. (2), 
44:412-431, 1938. 

[9] H. Davenport. On the product of three homogeneous linear forms. III. Proc. London Math. Soc. (2), 
45:98-125, 1939. 

[10] H. Davenport. Note on the product of three homogeneous linear forms. J. London Math. Soc, 16:98- 
101, 1941. 

[11] H. Davenport. On the product of three homogeneous linear forms. IV. Proc. Cambridge Philos. Soc, 
39:1-21, 1943. 

[12] J. W. Demmel. Applied numerical linear algebra. Society for Industrial and Applied Mathematics 

(SIAM), Philadelphia, PA, 1997. 
[13] K. Hessenberg. Thesis. Darmstadt, Germany: Technische Hochschule, 1942. 

[14] O. Karpenkov. On the triangulations of tori associated with two-dimensional continued fractions of 
cubic irrationalities. Fund. Anal. AppL, 38(2):102-110, 2004. Russian version: Funkt. Anal. Prilozh. 
38 (2), 2004, 28-37. 

[15] O. Karpenkov. On two-dimensional continued fractions of hyperbolic integer matrices with small 

norm. Russian Math. Surveys, 59(5):959-960, 2004. 
[16] O. Karpenkov. Constructing multidimensional periodic continued fractions in the sense of Klein. 

Math. Comp., 78(267):1687-1711, 2009. 
[17] O. Karpenkov. On determination of periods of geometric continued fractions for two-dimensional 

algebraic hyperbolic operators. Math. Notes, 88(1-2) :28-38, 2010. Russian version: Mat. Zametki, 

88(1), (2010), 30-42. 

[18] O. Karpenkov. Multidimensional Gauss Reduction Theory for conjugacy classes of SL(n,Z). Preprint, 
2012. 

[19] F. Klein. Uber eine geometrische Auffassung der gewohnliche Kettenbruchentwicklung. Nachr. Ges. 

Wiss. Gdttingen, Math-phys. Kl, 3:352-357, 1895. 
[20] F. Klein. Sur une representation gcometrique de developpement en fraction continue ordinaire. Nouv. 

Ann. Math., 15(3):327-331, 1896. 
[21] M. L. Kontsevich and Yu. M. Suhov. Statistics of Klein polyhedra and multidimensional continued 

fractions. In Pseudoperiodic topology, volume 197 of Amer. Math. Soc. Transl. Ser. 2, pages 9-27. 

Amer. Math. Soc, Providence, RI, 1999. 
[22] E. I. Korkina. Two-dimensional continued fractions. The simplest examples. Trudy Mat. Inst. 

Steklov., 209(Osob. Gladkikh Otobrazh. s Dop. Strukt.):143-166, 1995. 
[23] E. I. Korkina. The simplest 2-dimcnsional continued fraction. J. Math. Sci., 82(5):3680-3685, 1996. 

Topology, 3. 

[24] G. Lachaud. Polyedre d'Arnol'd et voile d'un cone simplicial: analogues du theoreme de Lagrange. 

C. R. Acad. Sci. Paris Ser. I Math., 317(8):711-716, 1993. 
[25] G. Lachaud. Voiles et polyhedres de Klein. Act. Sci. Ind., Hermann, 2002. 

[26] J. Lewis and D. Zagier. Period functions and the Selberg zeta function for the modular group. In The 
mathematical beauty of physics (Saclay, 1996), volume 24 of Adv. Ser. Math. Phys., pages 83-97. 
World Sci. Publ., River Edge, NJ, 1997. 



24 



OLEG KARPENKOV 



[27] Y. I. Manin and M. Marcolli. Continued fractions, modular symbols, and noncommutative geometry. 

Selecta Math. (N.S.), 8(3):475-521, 2002. 
[28] A. Markoff. Sur les formes quadratiques binaires indefinies. Math. Ann., 15(3-4) :381-406, 1879. 
[29] J.-O. Moussafir. Voiles et Polyedres de Klein: Geometrie, Algorithmes et Statistiques. docteur en 

sciences these, Universite Paris IX - Dauphine, 2000. 
[30] J. M. Ortega and H. F. Kaiser. The LL T and QR methods for symmetric tridiagonal matrices. 

Comput. J., 6:99-101, 1963/1964. 
[31] J. Stoer and R. Bulirsch. Introduction to numerical analysis, volume 12 of Texts in Applied Mathe- 
matics. Springer- Verlag, New York, third edition, 2002. Translated from the German by R. Bartels, 

W. Gautschi and C. Witzgall. 
[32] L. N. Trefethen and D. Bau. Numerical linear algebra. Society for Industrial and Applied Mathematics 

(SIAM), Philadelphia, PA, 1997. 
[33] G. F. Voronoi. On a Generalization of the Algorithm of Continued Fraction. Collected works in three 

volumes (In Russian). USSR Ac. Sci., Kiev., 1952. 

E-mail address, Oleg Karpcnkov: karpenk@mccme.ru 

TU Graz /Kopernikusgasse 24, A 8010 Graz, Austria/ 



